String-net condensation: A physical mechanism for topological phases 
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We show that quantum systems of extended objects naturally give rise to a large class of ex- 
otic phases - namely topological phases. These phases occur when the extended objects, called 
"string- nets" , become highly fluctuating and condense. We derive exactly soluble Hamiltonians 
for 2D local bosonic models whose ground states are string-net condensed states. Those ground 
states correspond to 2D parity invariant topological phases. These models reveal the mathemat- 
ical framework underlying topological phases: tensor category theory. One of the Hamiltonians - 
a spin- 1/2 system on the honeycomb lattice - is a simple theoretical realization of a fault tolerant 
quantum computer. The higher dimensional case also yields an interesting result: we find that 3D 
string-net condensation naturally gives rise to both emergent gauge bosons and emergent fermions. 
Thus, string-net condensation provides a mechanism for unifying gauge bosons and fermions in 3 
and higher dimensions. 

PACS numbers: 11.15.-q, 71.10.-w 



I. INTRODUCTION 

For many years, it was thought that Landau's theory 
of symmetry breaking [1] could describe essentially all 
phases and phase transitions. It appeared that all con- 
tinuous phase transitions were associated with a broken 
symmetry. However, after the discovery of the fractional 
quantum Hall (FQH) effect, it was realized that FQH 
states contain a new type of order - topological order 
- that is beyond the scope of Landau theory (for a re- 
view, see Ref. [2]). Since then the study of topological 
phases in condensed matter systems has been an active 
area of research. Topological phases have been investi- 
gated in a variety of theoretical and experimental sys- 
tems, ranging from FQH systems [3~6], quantum dimer 
models [7-10] , quantum spin models [11-19], to quan- 
tum computing [20, 21], or even superconducting states 
[22, 23]. This work has revealed a host of interesting 
theoretical phenomena and applications, including frac- 
tionalization, anyonic quasiparticles, and fault tolerant 
quantum computation. Yet, a general theory of topolog- 
ical phases is lacking. 

One way to reveal the gaps in our understanding is 
to compare with Landau's theory of symmetry break- 
ing phases. Landau theory is based on (a) the physical 
concepts of long range order, symmetry breaking, and or- 
der parameters, and (b) the mathematical framework of 
group theory. These tools allow us to solve three impor- 
tant problems in the study of ordered phases. First, they 
provide low energy effective theories for general ordered 
phases: Ginzburg-Landau field theories [24]. Second, 
they lead to a classification of symmetry-breaking states. 
For example, we know that there are only 230 different 
crystal phases in three dimensions. Finally, they allow us 
to determine the universal properties of the quasiparticle 
excitations (e.g. whether they are gapped or gapless). In 
addition. Landau theory provides a physical picture for 
the emergence of ordered phases - namely particle con- 
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densation. 

Several components of Landau theory have been suc- 
cessfully reproduced in the theory of topological phases. 
For example, the low energy behavior of topological 
phases is relatively well understood on a formal level: 
topological phases are gapped and are described by topo- 
logical quantum field theories (TQFT's).[25] The prob- 
lem of physically characterizing topological phases has 
also been addressed. Ref. [2] investigated the "topologi- 
cal order" (analogous to long range order) that occurs 
in topological phases. The author showed that topo- 
logical order is characterized by robust ground state de- 
generacy, nontrivial particle statistics, and gapless edge 
excitations. [3, 13, 26] These properties can be used to 
partially classify topological phases. Finally, the quasi- 
particle excitations of topological phases have been ana- 
lyzed in particular cases. Unlike the symmetry breaking 
case, the emergent particles in topologically ordered (or 
more generally, quantum ordered) states include (decon- 
fined) gauge bosons [27, 28] as well as fermions (in three 
dimensions) [29, 30] or anyons (in two dimensions) [31]. 
Fermions and anyons can emerge as collective excitations 
of purely bosonic models. 

Yet, the theory of topological phases is still incomplete. 
The theory lacks two important components: a physical 
picture (analogous to particle condensation) that clarifies 
how topological phases emerge from microscopic degrees 
of freedom, and a mathematical framework (analogous 
to group theory) for characterizing and classifying these 
phases. 

In this paper, we address these two issues for a large 
class of topological phases which we call "doubled" topo- 
logical phases. On a formal level, "doubled" topological 
phases are phases that are described by a sum of two 
TQFT's with opposite chiralities. Physically, they are 
characterized by parity and time reversal invariance. Ex- 
amples include all discrete lattice gauge theories, and all 
doubled Chern-Simons theories. It is unclear to what ex- 
tent our results generalize to chiral topological phases - 
such as in the FQH effect. 

We first address the problem of the physical picture 
for doubled topological phases. We argue that in these 
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t/U « 1 t/U » 1 

FIG. 1: A schematic phase diagram for the generic string-net 
Hamiltonian (3). When t/U (the ratio of the kinetic energy 
to the string tension) is small the system is in the normal 
phase. The ground state is essentially the vacuum with a few 
small string-nets. When t/U is largo the string-nets condense 
and large fluctuating string-nets fill all of space. We expect a 
phase transition between the two states at some t/U of order 
unity. We have omitted string labels and orientations for the 
sake of clarity. 

phases, local energetic constraints cause the microscopic 
degrees of freedom to organize into effective extended 
objects called "string-nets". At low energies, the mi- 
croscopic Hamiltonian effectively describes the dynamics 
of these extended objects. If the kinetic energy of the 
string-nets dominates the string-net tension, the string- 
nets "condense" : large string-nets with a typical size on 
the same order as the system size fill all of space (see 
Fig. 1). The result is a doubled topological phase. Thus, 
just as traditional ordered phases arise via particle con- 
densation, topological phases originate from "string-net 
condensation." 

This physical picture naturally leads to a solution 
to the second problem - that of finding a mathemati- 
cal framework for classifying and characterizing doubled 
topological phases. We show that each topological phase 
is associated with a mathematical object known as a 
"tensor category." [32] Here, we think of a tensor cat- 
egory as a 6 index object F^i'^^^ which satisfies certain 
algebraic equations (8). The mathematical object i^/^^ 
characterizes different topological phases and determines 
the universal properties of the quasiparticle excitations 
(e.g. statistics) just as the symmetry group does in Lan- 
dau theory. We feel that the mathematical framework 
of tensor categories, together with the physical picture 
of string-net condensation provides a general theory of 
(doubled) topological phases. 

Our approach has the additional advantage of provid- 
ing exactly soluble Hamiltonians and ground state wave 
functions for each of these phases. Those exactly solu- 
ble Hamiltonians describe local bosonic models (or spin 
models). They realize all discrete gauge theories (in any 
dimension) and all doubled Chern-Simons theories (in 
(2 -|- 1) dimensions). One of the Hamiltonians - a spin- 
1/2 model on the honeycomb lattice - is a simple theo- 
retical realization of a fault tolerant quantum computer 
[33]. The higher dimensional models also yield an in- 
teresting result: we find that (3 -|- 1)D string- net con- 



densation naturally gives rise to both emerging gauge 
bosons and emerging fcrmions. Thus, string-net conden- 
sation provides a mechanism for unifying gauge bosons 
and fermions in (3 + 1) and higher dimensions. 

We feel that this constructive approach is one of the 
most important features of this paper. Indeed, in the 
mathematical community it is well known that topolog- 
ical field theory, tensor category theory and knot theory 
are all intimately related [34-36]. Thus it is not sur- 
prising that topological phases are closely connected to 
tensor categories and string-nets. The contribution of 
this paper is our demonstration that these elegant mathe- 
matical relations have a concrete realization in condensed 
matter systems. 

The paper is organized as follows. In sections II and 
HI, we introduce the string-net picture, first in the case 
of deconfined gauge theories, and then in the general 
case. We argue that all doubled topological phases are 
described by string-net condensation. 

The rest of the paper is devoted to developing a the- 
ory of string-net condensation. In section IV, we con- 
sider the case of (2-1-1) dimensions. In parts A and B, 
we construct string-net wave functions and Hamiltonians 
for each (2 -|- 1)D string-net condensed phase. Then, in 
part C, we use this mathematical framework to calculate 
the universal properties of the quasiparticle excitations 
in each phase. In section V, we discuss the generaliza- 
tion to 3 and higher dimensions. In the last section, we 
present several examples of string-net condensed states - 
including a spin- 1/2 model theoretically capable of fault 
tolerant quantum computation. The main mathematical 
calculations can be found in the appendix. 

II. STRING-NETS AND GAUGE THEORIES 

In this section, we introduce the string-net picture in 
the context of gauge theory [27, 37, 38]. We point out 
that all deconfined gauge theories can be understood as 
string-net condensates where the strings are essentially 
electric flux lines. We hope that this result provides in- 
tuition for (and motivates) the string-net picture in the 
general case. 

We begin with the simplest gauge theory - Z2 lattice 
gauge theory [39] . The Hamiltonian is 

i p edges of p 

where cr^'^'^ are the Pauli matrices, and I, i, p label the 
sites, links, and plaquettes of the lattice. The Hilbert 
space is formed by states satisfying 

n ^?i*)=i^)' (2) 

legs of / 

for every site /. For simplicity we will restrict our dis- 
cussion to trivalent lattices such as the honeycomb lattice 
(see Fig. 2). 

It is well known that Z2 lattice gauge theory is dual 
to the Ising model in (2 -|- 1) dimensions [40]. What 
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FIG. 2: The constraint term Yli^ s oC i '^i ^^'^ magnetic term 
ricdgcs of p '^j -^2 lattice gauge theory. In the dual picture, 
we regard the links with = —1 as being occupied by a 
string, and the links with ct"^ = +1 as being unoccupied. The 
constraint term then requires the strings to be closed - as 
shown on the right. 

is less well known is that there is a more general dual 
description of Z2 gauge theory that exists in any number 
of dimensions [41]. To obtain this dual picture, we view 
links with cr^ = — 1 as being occupied by a string and 
links with cr^ — +1 as being unoccupied. The constraint 
(2) then implies that only closed strings are allowed in 
the Hilbert space (Fig. 2). 

In this way, Z2 gauge theory can be reformulated as a 
closed string theory, and the Hamiltonian can be viewed 
as a closed string Hamiltonian. The electric and mag- 
netic energy terms have a simple interpretation in this 
dual picture: the "electric energy" —UJ2i'^i is a string 
tension while the "magnetic energy" t J2p Ilodgcs of p 
is a string kinetic energy. The physical picture for the 
confining and deconfined phases is also clear. The con- 
fining phase corresponds to a large electric energy and 
hence a large string tension U ^ t. The ground state 
is therefore the vacuum configuration with a few small 
strings. The deconfined phase corresponds to a large 
magnetic energy and hence a large kinetic energy. The 
ground state is thus a superposition of many large string 
configurations. In other words, the deconfined phase of 
Z2 gauge theory is a quantum liquid of large strings - a 
"string condensate." (Fig. 3a). 

A similar, but more complicated, picture exists for 
other deconfined gauge theories. The next layer of com- 
plexity is revealed when we consider other Abelian the- 
ories, such as C/(l) gauge theory. As in the case of Z2, 
U{1) lattice gauge theory can be reformulated as a the- 
ory of electric flux lines. However, unlike Z2, there is 
more then one type of flux line. The electric flux on a 
link can take any integral value in U{1) lattice gauge the- 
ory. Therefore, the electric flux lines need to be labeled 
with integers to indicate the amount of flux carried by 
the line. In addition, the flux lines need to be oriented 
to indicate the direction of the flux. The final point is 
that the flux lines don't necessarily form closed loops. 
It is possible for three flux lines Ei, E2, to meet at a 
point, as long as Gauss' law is obeyed: E1+E2+E3 — 0. 
Thus, the dual formulation of U{1) gauge theory involves 
not strings, but more general objects: networks of strings 
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FIG. 3: Typical string-net configurations in the dual formu- 
lation of (a) Z2, (b) (7(1), and (c) SU{2) gauge theory. In 
the case of (a) Z2 gauge theory, the string-net configura- 
tions consist of closed (non-intersecting) loops. In (b) U{1) 
gauge theory, the string-nets are oriented graphs with edges 
labeled by integers. The string-nets obey the branching rules 
El + E2 + E3 = for any three edges meeting at a point. 
In the case of (c) SU{2) gauge theory, the string-nets con- 
sist of (unoriented) graphs with edges labeled by half-integers 
1/2, 1, 3/2, .... The branching rules are given by the triangle 
inequality: {Ei,E2, E3} are allowed to meet at a point if and 
only if El < E2 + E3, E2 < E3 + Ei, E3 < Ei + E2, and 
El + E2 + E3 is an integer. 



(or "string-nets" ) . The strings in a string- net are labeled, 
oriented, and obey branching rules, given by Gauss' law 
(Fig. 3b). 

This "string-net" picture exists for general gauge the- 
ories. In the general case, the strings (electric flux lines) 
are labeled by representations of the gauge group. The 
branching rules (Gauss' law) require that if three strings 
El, E2, i?3 meet at a point, then the product of the rep- 
resentations El ® E2® E3 must contain the trivial rep- 
resentation. (For example, in the case of SU{2), the 
strings are labeled by half- integers E — 1/2,1,3/2,..., 
and the branching rules are given by the triangle inequal- 
ity: {El, E2, E3} are allowed to meet at a point if and 
only if El < E2 + E3, E2 < E3 + Ei, E3 < Ei + E2 and 
E1+E2 + E'i is an integer (Fig. 3c)) [37] . These string- 
nets provide a general dual formulation of gauge theory. 
As in the case of Z2, the deconfined phase of the gauge 
theory always corresponds to highly fluctuating string- 
nets - a string-net condensate. 



III. GENERAL STRING-NET PICTURE 

Given the large scope of gauge theory, it is natural 
to wonder if string-nets can describe more general topo- 
logical phases. In this section we will discuss this more 
general string-net picture. (Actually, we will not discuss 
the most general string-net picture. We will focus on a 
special case for the sake of simplicity. See Appendix A 
for a discussion of the most general picture). 

We begin with a more detailed deflnition of "string- 
nets." As the name suggests, string-nets are networks of 
strings. We will focus on trivalent networks where each 
node or branch point is attached to exactly 3 strings. The 
strings in a string-net are oriented and come in various 
"types." Only certain combinations of string types are 
allowed to meet at a node or branch point. To specify 
a particular string-net model, one needs to provide the 
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FIG. 4: The orientation convention for the branching rules. 



FIG. 5: i and i* label strings with opposite orientations. 

following data: 

1. String types: The number of different string types 
N. For simplicity, we will label the different string 
types with the integers i = 1, ...,N. 

2. Branching rules: The set of all triplets of string- 
types {{i,j,k}...} that are allowed to meet at a 
point. (See Fig. 4). 

3. String orientations: The dual string type i* as- 
sociated with each string type i. The duality must 
satisfy (i*)* = i. The type-z* string corresponds to 
the type-i string with the opposite orientation. If 
i = i*, then the string is unoriented (See Fig. 5). 

This data describes the detailed structure of the string- 
nets. The Hilbert space of the string-net model is then 
defined in the natural way. The states in the Hilbert 
space are simply linear superpositions of different spatial 
configurations of string-nets. 

Once the Hilbert space has been specified, we can 
imagine writing down a string-net Hamiltonian. The 
string- net Hamiltonian can be any local operator which 
acts on quantum string- net states. A typical Hamiltonian 
is a sum of potential and kinetic energy pieces: 

H = UHu + tHt (3) 

The kinetic energy Ht gives dynamics to the string-nets, 
while the potential energy Hjj is typically some kind of 
string tension. When U >> t, the string tension domi- 
nates and we expect the ground state to be the vacuum 
state with a few small string-nets. On the other hand, 
when t >> U , the kinetic energy dominates, and we 
expect the ground state to consist of many large fluc- 
tuating string- nets. We expect that there is a quantum 
phase transition between the two states at some t/U on 
the order of unity. (See Fig. 1). Because of the analogy 
with particle condensation, we say that the large t, highly 
fluctuating string-net phase is "string-net condensed." 

This notion of string-net condensation provides a nat- 
ural physical mechanism for the emergence of topological 
phases in real condensed matter systems. Local energetic 
constraints can cause the microscopic degrees of freedom 
to organize into effective extended objects or string-nets. 
If the kinetic energy of these string-nets is large, then 
they can condense giving rise to a topological phase. The 
type of topological phase is determined by the structure 
of the string-nets, and the form of string-net condensa- 
tion. 



But how general is this picture? In the previous sec- 
tion, we pointed out that all dcconfincd gauge theories 
can be viewed as string-net condensates. In fact, mathe- 
matical results suggest that the string-net picture is even 
more general. In (2 + 1) dimensions, all so-called "dou- 
bled" topological phases can be described by string-net 
condensation (provided that we generalize the string-net 
picture as in Appendix A). [34] Physically, this means 
that the string-net picture can be applied to essentially 
all parity and time reversal invariant topological phases 
in (2-1-1) dimensions. Examples include all discrete gauge 
theories, and all doubled Chern-Simons theories. The sit- 
uation for dimension c? > 2 is less well understood. How- 
ever, we know that string-net condensation quite gener- 
ally describes all lattice gauge theories with or without 
emergent Fermi statistics. 



IV. STRING-NET CONDENSATION IN (2 + 1) 
DIMENSIONS 

A. Fixed-point wave functions 

In this section, we attempt to capture the imiversal 
features of string- net condensed phases in (2 -|- 1) dimen- 
sions. Our approach, inspired by Ref. [35-37, 42-44], 
is based on the string-net wave function. We construct 
a special "fixed-point" wave function for each string-net 
condensed phase. We believe that these "fixed-point" 
wave functions capture the universal properties of the 
corresponding phases. Each "fixed-point" wave function 
is associated with a six index object F^^^ that satisfies 
certain algebraic equations (8). In this way, we derive 
a one-to-one correspondence between doubled topologi- 
cal phases and tensor categories F^^^. We would like 
to mention that a related result on the classification of 
{2 + l)D topological quantum field theories was obtained 
independently in the mathematical community. [34] 

Let us try to visualize the wave function of a string-net 
condensed state. Though we haven't defined string-net 
condensation rigorously, we expect that a string-net con- 
densed state is a superposition of many different string- 
net configurations. Each string-net configuration has 
a size typically on the same order as the system size. 
The large size of the string-nets implies that a string-net 
condensed wave function has a non-trivial long distance 
structure. It is this long distance structure that distin- 
guishes the condensed state from the "normal" state. 

In general, we expect that the universal features of a 
string-net condensed phase are contained in the long dis- 
tance character of the wave functions. Imagine compar- 
ing two different string-net condensed states that belong 
to the same quantum phase. The two states will have 
different wave functions. However, by the standard RG 
reasoning, we expect that the two wave functions will 
look the same at long distances. That is, the two wave 
functions will only differ in short distance details - like 
those shown in Fig. 6. 

Continuing with this line of thought, we imagine per- 
forming an RG analysis on ground state functions. All 
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FIG. 6: Three pairs of string-net configurations that difi^er 

only in their short distance structure. We expect string-not 
wave functions in the same quantum phase to only differ by 
these short distance details. 




FIG. 7: A schematic RG flow diagram for a string-net model 
with 4 string-net condensed phases a. b, c, and d. All the 
states in each phase flow to fixed-points in the long distance 
limit. The corresponding fixed-point wave functions $a, $6, 
3>c, and $d capture the universal long distance features of 
the associated quantum phases. Our ansatz is that the fixed- 
point wave functions "1> are described by local constraints of 
the form (4-7). 

the states in a string-net condensed phase should flow 
to some special "fixed-point" state. We expect that the 
wave function of this state captures the universal long 
distance features of the whole quantum phase. (See Fig. 

In the following, we will construct these special fixed- 
point wave functions. Suppose $ is some fixed-point 
wave function. We know that $ is the ground state of 
some fixed-point Hamiltonian H. Based on our experi- 
ence with gauge theories, we expect that H is free. That 
is, if is a sum of local string kinetic energy terms with 
no string tension terms: 

H = tHt = t Ht^i 

i 

In particular, H is unfrustrated, and the ground state 
wave function minimizes the expectation values of all 
the kinetic energy terms {Ha} simultaneously. Minimiz- 
ing the expectation value of an individual kinetic energy 
term Ht i is equivalent to imposing a local constraint 
on the ground state wave function, namely = 
(where Ei is the smallest eigenvalue of Hij). 
We conclude that the wave function $ can be specified 
uniquely by local constraint equations. The local con- 
straints are linear relations between several string-net 
amplitudes ^{Xi),^{X2),^{Xs)--- where the configura- 
tions Xi,X2,X3... only differ by local transformations. 

To derive these local constraints from first principles 
is difficult, so we will use a more heuristic approach. We 
will first guess the form of the local constraints {ie guess 
the form of the fixed-point wave function). Then, in the 



next section, we will construct the fixed-point Hamilto- 
nian and show that its ground state wave function does 
indeed satisfy these local relations. Our ansatz is that 
the local constraints can be put in the following graphi- 
cal form: 

$ ^1 O'^ =d^<^ (^1^ (5) 
* (Bi^^jB) ='^^.-* (~tOT ) (6) 

* {mHm) = E * {mSJm) 

Here, i, j, k etc. are arbitrary string types and the 
shaded regions represent arbitrary string-net c;onfigura- 
tions. The di are complex numbers. The 6 index sym- 
bol F^j^ is a complex numerical constant that depends 
on 6 string types i, j, m, fc, I, and n. If one or more 
of the branchings {i,j, to}, {k, I, m*}, {i, n, I}, {j, k, n*} is 
illegal, the value of the symbol F^j^ is unphysical. How- 
ever, for simplicity, we will set = in this case. 

The local rules (4-7) are written using a new nota- 
tional convention. According to this convention, the in- 
dices i,j,k etc., can take on the value i = in addition 
to the N physical string types z = 1, ...N. We think of 
the i = string as the "empty string" or "null string." 
It represents empty space - the vacuum. Thus, we can 
convert labeled string-nets to our old convention by sim- 
ply erasing all the z = strings. The branching rules and 
dualities associated with z = are defined in the obvious 
way: 0* = 0, and {i, j, 0} is allowed if and only if z = j* . 
Our convention serves two purposes: it simplifies nota- 
tion (each equation in (4- 7) represents several equations 
with the old convention) , and it reveals the mathematical 
framework underlying string-net condensation. 

We now briefly motivate these rules. The flrst rule (4) 
constrains the wave function $ to be topologically invari- 
ant. It requires the quantum mechanical amplitude for 
a string-net configuration to only depend on the topol- 
ogy of the configuration: two configurations that can be 
continuously deformed into one another must have the 
same amplitude. The motivation for this constraint is 
our expectation that topological string-net phases have 
topologically invariant fixed-points. 

The second rule (5) is motivated by the fundamental 
property of RG fixed-points: scale invariancc. The wave 
function $ should look the same at all distance scales. 
Since a closed string disappears at length scales larger 
then the string size, the amplitude of an arbitrary string- 
net configuration with a closed string should be propor- 
tional to the amplitude of the string-net configuration 
alone. 

The third rule (6) is similar. Since a "bubble" is irrel- 
evant at long length scales, we expect 
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But if i ^ J, the configuration ' ' is not allowed 
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0. Wc conclude that the amplitude 



for the bubble configuration vanishes when i ^ j (6). 

The last rule is less well-motivated. The main point 
is that the first three rules are not complete: another 
constraint is needed to specify the ground state wave 
function uniquely. The last rule (7) is the simplest lo- 
cal constraint with this property. An alternative motiva- 
tion for this rule is the fusion algebra in conformal field 
theory. [45] 

The local rules (4-7) uniquely specify the fixed-point 
wave function i>. The universal features of the string- net 

condensed state arc captured by these rules. Equiva- 
lently, they are captured by the six index object F^j™ , 
and the numbers di. 

However, not every choice of (-F^;™, di) corresponds to 
a string-net condensed phase. In fact, a generic choice 
of (-Pfcf™, di) will lead to constraints (4-7) that are not 
self-consistent. The only {F^j^ , di) that give rise to self- 
consistent rules and a well-defined wave function $ are 
(up to a trivial rescaling) those that satisfy 



ijk _ 



JV 

E 

n=0 



pijm _ plkm* 
kin jin 



pmlq pjip pjs'i 
kp'n mns'^lkr* 



ViVj 



Oijk 



^Ikn* 



pimj '^mUn 
*^*"' VjVl 



pjip pnq 
q*kr* mis* 



(8) 



where Vi = Vi» = \fdi (and = 1). (See appendix B). 
Here, we have introduced a new object 5ijk defined by 
the branching rules: 



if k} is allowed, 
otherwise. 



(9) 



There is a one-to-one correspondence between (2-\-l)D 
string- net condensed phases and solutions of (8). These 
sohitions correspond to mathematical objects known as 
tensor categories. [32] Tensor category theory is the fun- 
damental mathematical framework for string-net conden- 
sation, just as group theory is for particle condensation. 
We have just shown that it gives a complete classifica- 
tion of (2 -|- 1)D string- net condensed phases (or equiva- 
lently doubled topological phases): each phase is associ- 
ated with a different solution to (8). We will show later 
that it also provides a convenient framework for deriving 
the physical properties of quasiparticles. 

It is highly non- trivial to find solutions of (8). How- 
ever, it turns out each group G provides a solution. The 
solution is obtained by (a) letting the string-type index i 
run over the irreducible representations of the group, (b) 
letting the numbers di be the dimensions of the represen- 
tations and (c) letting the 6 index object -F^;™ be the 6j 
symbol of the group. The low energy effective theory of 
the corresponding string-net condensed state turns out to 
be a deconfined gauge theory with gauge group G. An- 
other class of solutions can be obtained from &j symbols 



of quantum groups. It turns out that in these cases, the 
low energy effective theories of the corresponding string- 
net condensed states are doubled Chern-Simons gauge 
theories. These two classes of solutions are not necessar- 
ily exhaustive: Eq. (8) may have solutions other then 
gauge theories or Chern-Simons theories. Nevertheless, 
it is clear that gauge bosons and gauge groups emerge 
from string-net condensation in a very natural way. 

In fact, string- net condensation provides a new per- 
spective on gauge theory. Traditionally, we think of 
gauge theories geometrically. The gauge field is anal- 
ogous to an affine connection, and the field strength F^^i, 
is essentially a curvature tensor. From this point of view, 
gauge theory describes the dynamics of certain geometric 
objects (e.g. fiber bimdles). The gauge group determines 
the structure of these objects and is introduced by hand 
as part of the basic definition of the theory. In con- 
trast, according to the string- net condensation picture, 
the geometrical character of gauge theory is not funda- 
mental. Gauge theories are fundamentally theories of 
extended objects. The gauge group and the geometrical 
gauge structure emerge dynamically at low energies and 
long distances. A string- net system "chooses" a particu- 
lar gauge group, depending on the coupling constants in 
the underlying Hamiltonian: these parameters determine 
a string-net condensed phase which in turn determines a 
solution to (8). The nature of this solution determines 
the gauge group. 

One advantage of this alternative picture is that it uni- 
fies two seemingly unrelated phenomena: gauge interac- 
tions and Fermi statistics. Indeed, as we will show in 
section V, string-net condensation naturally gives rise to 
both gauge interactions and Fermi statistics (or fractional 
statistics in (2 + 1)D). In addition, these structures al- 
ways appear together. [29] 



B. Fixed-point Hamiltonians 

In this section, we construct exactly soluble lattice spin 
Hamiltonians with the fixed-point wave functions $ as 
ground states. These Hamiltonians provide an explicit 
realization of all (2 -|- 1)D string-net condensates and 
therefore all (2 + 1)D doubled topological phases (pro- 
vided that we generalize these models as discussed in 
Appendix A). In the next section, we will use them to 
calculate the physical properties of the quasiparticle ex- 
citations. 

For every (F^;™, rfi) satisfying the self-consistency con- 
ditions (8) and the unitarity condition (14), we can con- 
struct an exactly soluble Hamiltonian. Let us first de- 
scribe the Hilbert space of the exactly soluble model. The 
model is a spin system on a (2D) honeycomb lattice, with 
a spin located on each link of the lattice. Each "spin" can 
be in + 1 different states labeled by i = 0, 1, A^. We 
assign each link an arbitrary orientation. When a spin 
is in state i, we think of the link as being occupied by 
a type-i string oriented in the appropriate direction. We 
think of the type-0 string or null string as the vacuum 
{ie no string on the link). 
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FIG. 8: A picture of the lattice spin model (10). The electric 
charge operator Qi acts on the three spins adjacent to the 
vertex /, while the magnetic energy operator Bp acts on the 
12 spins adjacent to the hexagonal plaquette p. The term Qi 
constrains the string-nets to obey the branching rules, while 
Bp provides dynamics. A typical state satisfying the low- 



energy constraints is shown on the right, 
have spins in the i = state. 



The empty links 



Bt 




Then the operators are defined by 




B 



(abcdef) 



s,g h i j k I 
p,ghijkl 

pal* g pbg'h pen i pu-i j pcj 



s* j'i'* s* k' j'* s*l'k'' 



(12) 



(13) 



(See appendix C for a graphical representation of Bp). 
One can check that the Hamiltonian (10) is Hermitian if 
F satisfies 



pi J m 
^k'l'n* 



i^kln ) 



(14) 



by 



The exactly soluble Hamiltonian for our model is given 



N 



(10) 



where the sums run over vertices I and plaquettes p of 
the honeycomb lattice. The coefficients as satisfy Os* — 
a* but are otherwise arbitrary. 

Let us explain the terms in (10). We think of the 
first term Qz as an electric charge operator. It measures 
the "electric charge" at site I, and favors states with no 
charge. It acts on the 3 spins adjacent to the site /: 



Qi 



i j 



(11) 



where Sijk is the branching rule symbol (9). Clearly, this 
term constrains the strings to obey the branching rules 
described by Sij^.. With this constraint the low energy 
Hilbert space is essentially the set of all allowed string- 
net configurations on a honeycomb lattice. (See Fig. 8). 

We think of the second term Bp as a magnetic flux 
operator. It measures the "magnetic flux" though the 
plaquette p (or more precisely, the cosine of the magnetic 
flux) and favors states with no flux. This term provides 
dynamics for the string-net configurations. 

The magnetic flux operator Bp is a linear combination 
of (A^-l-l) terms Bp, s = 0, 1, N. Each Bp is an opera- 
tor that acts on the 12 links that are adjacent to vertices 
of the hexagon p. (See Fig. 8). Thus, the Bp are essen- 
tially (7V-|-1)^^ X (TV -1-1)^^ matrices. However, the action 
of Bp does not change the spin states on the 6 outer links 
of p. Therefore the Bp can be block diagonalized into 
{N + If blocks, each of dimension {N + 1)^ x {N + if. 

= 0, l,...,iV, de- 
1)^ matrices: 



^pHhijki ^ ' (abcdef ), with a, 6, c. 
note the matrix elements of these {N 



in addition to (8). Our model is only applicable to topo- 
logical phases satisfying this additional constraint. We 
believe that this is true much more generally: only topo- 
logical phases satisfying the unitarity condition (14) are 
physically realizable. 

The Hamiltonian (10) has a number of interesting 
properties, provided that {Fljl^ , di) satisfy the self- 
consistency conditions (8). It turns out that: 



1. The Bp and Q/'s all commute with each other. 
Thus the Hamiltonian (10) is exactly soluble. 

2. Depending on the choice of the coefficients Os, the 
system can be in iV -I- 1 different quantum phases. 

3. The choice ag = ^ corresponds to a topolog- 

ical phase with a smooth continuum limit. The 
ground state wave function for this parameter 
choice is topologically invariant, and obeys the lo- 
cal rules (4-7). It is precisely the wave function $, 
defined on a honeycomb lattice. Furthermore, Q/, 
Bp are projection operators in this case. Thus, the 
ground state satisfies Qi = Bp = 1 for all /, p, 
while the excited states violate these constraints. 

The Hamiltonian (10) with the above choice of pro- 
vides an exactly soluble realization of the doubled topo- 
logical phase described by Fjfi^ . We can obtain some in- 
tuition for the Hamiltonian (10) by considering the case 
where -F^;™ is the 6j symbol of some group G. In this 
case, it turns out that Qi and Bp are precisely the elec- 
tric charge and magnetic fiux operators in the standard 
lattice gauge theory with group G. Thus, (10) is the 
usual Hamiltonian of lattice gauge theory, except with 
no electric field term. This is nothing more than the 
well-known exactly soluble Hamiltonian of lattice gauge 
theory. [20, 39] In this way, our construction can be 
viewed as a natural generalization of lattice gauge the- 
ory. 
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In this paper, we will focus on the smooth topolog- 
ical phase corresponding to the parameter choice = 
^2 (see appendix C). However, we would like to 

mention that the other N quantum phases also have non- 
trivial topological (or quantum) order. However, in these 
phases, the ground state wave function does not have a 
smooth continuum limit. Thus, these arc new topological 
phases beyond those described by continuum theories. 



trix element of a general typc-s simple string operator 
W{P) between an initial spin state ii, ...ipf and final spin 
state i[, ...i'j^ is of the form 



N 



W. 



(15) 



Vfe=i 



\k=l 



where ei, cat are the spin states of the N "legs" of P 
(see Fig. 9) and 



C. Quasiparticle excitations 

In this section, wc find the quasiparticle excitations 
of the string-net Hamiltonian (10), and calculate their 
statistics (e.g. the twists 6a and the S matrix Sa/?)- We 

will (mly consider the topological phase with smooth con- 
tinuum limit. That is, we will choose ag = 2 in our 

lattice model. 

Recall that the ground state satisfies Qi = Bp = 1 for 
all vertices /, and all plaquettes p. The quasiparticle ex- 
citations correspond to violations of these constraints for 
some local collection of vertices and plaquettes. We are 
interested in the topological properties (e.g. statistics) of 
these excitations. 

We will focus on topologically nontrivial quasiparticles 
- that is, particles with nontrivial statistics or mutual 
statistics. By the analysis in Rcf. [29] , wc know that these 
types of particles are always created in pairs, and that 
their pair creation operator has a string-like structure, 
with the newly created particles appearing at the ends. 
(See Fig. 9). The position of this string operator is 
unobservable in the string-net condensed state - only the 
endpoints of the string arc observable. Thus the two ends 
of the string behave like independent particles. 

If the two endpoints of the string coincide so that the 
string forms a loop, then the associated closed string 
operator commutes with the Hamiltonian. This follows 
from the fact that the string is truly unobservable; the 
action of an open string operator on the ground state 
depends only on its endpoints. 

Thus, each topologically nontrivial quasiparticle is as- 
sociated with a (closed) string operator that commutes 
with the Hamiltonian. To find the quasiparticles, we need 
to find these closed string operators. 

An important class of string operators are what we 
will call "simple" string operators. The defining property 
of simple string operators is their action on the vacuum 
state. If wc apply a type-s simple string operator W{P) 
to the vacuum state, it creates a type-s string along the 
path of the string, P. Wc already have some examples of 
these operators, namely the magnetic fiux operators Bp. 
When Bp acts on the vacuum configuration |0), it creates 
a type-s string along the boundary of the plaquette p. 
Thus, we can think of Bp as a short type-s simple string 
operator, W{dp). 

We would like to construct simple string operators 
W{P) for arbitrary paths P = Ji,...,Jjv on the hon- 
eycomb lattice. Using the definition of Bp as a guide, we 
make the following ansatz. The string operator W{P) 
only changes the spin states along the path P. The ma- 



pek^yk-l ^ -J p ^^j^j^g ^^^^ 

F« = { (16) 
Kt' i'' ^^ if -P turns right at Jfe 



% if P turns right, left at Ik, Ik+i 



u)k - ■{ lliklii^'^k^ if p turns left, right at Ik, Jfe+i 

k 

1 , otherwise 

(17) 

Here, ujj, uj^j are two (complex) two index objects that 
characterize the string W. 

Note the similarity to the definition of Bp. The major 

difference is the additional factor OfcLi'-^fc- We conjec- 
ture that Iw^^^^^l = 1 for a type-s string, so Jl^i^fe 
is simply a phase factor that depends on the initial and 



final spin states ii, Z2, 



In, '-2j *JV 



This phase 



vanishes for paths P that make only left or only right 
turns, such as plaquette boundaries dp. In that case, the 
definition of W(P) coincides with Bp. 

A straightforward calculation shows that the opera- 
tor W{P) defined above commutes with the Hamiltonian 
(10) if i^pi^j satisfy 



3 Kjm* ■ 



iVjVs 



N 



Epji'k ^.npjl'n 
s*nl* k ksm* 

N 

E k pis'k 



(18) 



k=0 



The solutions to these equations give all the type-s simple 
string operators. 

For example, consider the case of Abelian gauge the- 
ory. In this case, the solutions to (18) can be divided 
into three classes. The first class is given by s 7^ 0, 
'^i ^IT^ ~ 'if^ ~ -i- These string operators create 
electric flux lines and the associated quasiparticles are 
electric charges. In more traditional nomenclature, these 
are known as (Wegner-) Wilson loop operators [39, 46]. 
The second class of solutions is given by s = 0, and 
ujj^^jj^ = {ujf^^^y* 1. These string operators cre- 
ate magnetic fiux lines and the associated quasiparti- 
cles are magnetic fluxes. The third class has s 7^ and 
= t^JlHv^Y u These strings create both elcc- 

trie and magnetic flux and the associated quasiparticles 
are electric charge/magnetic flux bound states. This ac- 
counts for all the quasiparticles in (2-|-l)_D Abelian gauge 
theory. Therefore, all the string operators are simple in 
this case. 
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However, this is not true for non-Abelian gauge the- 
ory or other (2 + l)D topological phases. To compute 
the quasiparticle spectrum of these more general theo- 
ries, we need to generalize the expression (15) for W{P) 
to include string operators that are not simple. 

One way to guess the more general expression for 
W{P) is to consider products of simple string operators. 
Clearly, if Wi{P) and W2{P) commute with the Hamil- 
tonian, then W{P) — Wi{P)-W2{P) also commutes with 
the Hamiltonian. Thus, we can obtain other string oper- 
ators by taking products of simple string operators. In 
general, the resulting operators are not simple. If Wi and 
W2 are type-si and type-S2 simple string operators, then 
the action of the product string on the vacuum state is: 



where \s) denotes the string state with a type-s string 
along the path P and the vacuum everywhere else. If we 
take products of more then two simple string operators 
then the action of the product string on the vacuum is 
of the form W^(P)|0) — J2s "-sk) where are some non- 
negative integers. 

We now generalize the expression for W{P) so that it 
includes arbitrary products of simple strings. Let W be 
a product of simple string operators, and let Us be the 
non-negative integers characterizing the action of W on 
the vacuum: VF(P)|0) = X^s'^sls). Then, one can show 
that the matrix elements of W{P) are always of the form 



N 



{Sk} 



N 



Tr 



\k=l 



\k=l 



where 



(19) 




FIG. 9: Open and closed string operators for the lattice 
spin model (10). Open string operators create quasiparti- 
cles at the two ends, as shown on the left. Closed string 
operators, as shown on the right, commute with the Hamil- 
tonian. The closed string operator W{P) only acts non- 
trivially on the spins along the path P — 7i,72... (thick 
line), but its action depends on the spin states on the 
legs (thin lines). The matrix element between an initial 

state 11,12,... and a final state i'i,i2,... is W^^'Av.. (£162...) = 



pie string and Wl^l^"' [eie2- 



Tr(- 



'j^ "^ o)/...) for a type-s sim- 
\ls,^.■■■) for a general string. 



and Vl satisfy 

N 



N 



^ ^ ^rsj-^kjm* ^ 



I '^J'^s 



s=0 



^t'nl'^^rtk^krm* 

n=0 
N 

Y,^sH'F:Q. (22) 



k=0 



•Id, 



if P tm-ns right, left at Jfc, Jfc_|_i 

if P turns left, right at J^, Ik+i 
otherwise 

(20) 

and are two 4 index objects that characterize 

the string operator W. For any quadruple of string types 
i,j,s,t, i^ltj: ^Itj) (complex) rectangular matrices 
of dimension rig x rit- Note that type-so simple string 
operators correspond to the special case where Ug = Ss^s- 



In this case, the matrices ^Itj^ 
numbers, and we can identify 



^l^j reduce to complex 



stj 



jf^sso^tso J 



•jSssoStso- 



(21) 



As we mentioned above, products of simple string op- 
erators are always of the form (19). In fact, we be- 
lieve that all string operators are of this form. Thus, 
we will use (19) as an ansatz for general string operators 
in (2 -I- 1)D topological phases. This ansatz is compli- 
cated algebraically, but like the definition of B^, it has a 
simple graphical interpretation (see appendix D). 

A straightforward calculation shows that the closed 
string W{P) commutes with the Hamiltonian (10) if H. 



The solutions {flmj^m) to these equations give all the 
different closed string operators Wm- However, not all 
of these solutions are really distinct. Notice that two 
solutions (riijCli), (Sl2,f^2) can be combined to form a 
new solution {fl' , fj'): 



sti 



(23) 



This is not surprising: the string operator W' corre- 
sponding to (il'jfj') is simply the sum of the two op- 
erators corresponding to (r2i,2, ^^1,2): W = W\ -I- W2- 

Given this additivity property, it is natural to con- 
sider the "irreducible" solutions (ri^, £7^,) that cannot be 
written as a sum of two other solutions. Only the "ir- 
reducible" string operators Wa create quasiparticle-pairs 
in the usual sense. Reducible string operators W create 
superpositions of different strings - which correspond to 
superpositions of different quasiparticles. [51] 

To analyze a topological phase, one only needs to find 
the irreducible solutions iyia^^a) to (22). The number 
M of such solutions is always finite. In general, each so- 
lution corresponds to an irreducible representation of an 
algebraic object. In the case of lattice gauge theory, there 
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is one solution for every irreducible representation of the 
quantum double D(G) of the gauge group G. Similarly, 
in the case of doubled Chern-Simons theories there is one 
solution for each irreducible representation of a doubled 
quantum group. 

The structure of these irreducible string operators Wa 
determines all the universal features of the topological 
phase. The number M of irreducible string operators is 
the number of different kinds of quasiparticles. The fu- 
sion rules WaWp = J2^=i ^Ixp^i determine how bound 
states of type-a and type-/3 quasiparticles can be viewed 
as a superposition of other types of quasiparticles. 

The topological properties of the quasiparticles are also 
easy to compute. As an example, we now derive two par- 
ticularly fundamental objects that characterize the spins 
and statistics of quasiparticles: the M twists 9 a and the 
M X M S-matrix, [13, 25, 35, 47]. 

The twists 6^ are defined to be statistical angles of 
the type-a quasiparticles. By the spin-statistics theo- 
rem they are closely connected to the quasiparticle spins 
Sq.: e*^° = e^'^*''°. We can calculate 9a by comparing the 
quantum mechanical amplitude for the following two pro- 
cesses. In the first process, we create a pair of quasipar- 
ticles a, a (from the ground state), exchange them, and 
then annihilate the pair. In the second process, we create 
and then annihilate the pair without any exchange. The 
ratio of the amplitudes for these two processes is precisely 

The amplitude for each process is given by the expec- 
tation value of the closed string operator Wa for a par- 
ticular path P: 



A2 = 



(24) 



(25) 



Here, |<I>) denotes the ground state of the Hamiltonian 
(10). 

Let (ilajfiajna) bc the irreducible solution corre- 
sponding to the string operator Wa- The above 
two amplitudes can be then be expressed in terms of 
(p,a,^a,na) (scc appendix D): 



S 

A2 = na,sds 



(26) 
(27) 



Combining these results, we find that the twists are given 
by 



Ai E.rf^Tr(^70^,,,.) 



A2 



(28) 



Just as the twists 9a are related to the spin and statis- 
tics of individual particle types a, the elements of the 
S-matrix, Sap describe the mutual statistics of two parti- 
cle types a, (3. Consider the following process: We create 
two pairs of quasiparticles a, a, /3, braid a around /3, 
and then annihilate the two pairs. The element Sap is 
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FIG. 10: A three dimensional trivalent lattice, obtained by 
splitting the sites of the cubic lattice. We replace each vertex 
of the cubic lattice with 4 other vertices as shown above. 



defined to be the quantum mechanical amplitude A of 
this process, divided by a proportionality factor D where 



The amplitude A can be cal- 



culated from the expectation value of Wa , W^ for two 
"linked" paths P: 



A 



(29) 



Expressing A in terms of (ila, SIq, n^), we find 

^o.p = ^ = ),Y. Tr(f^l.»,. ) • '^A'^'h.m' )d^dj (30) 



ijk 



V. STRING-NET CONDENSATION IN (3 + 1) 
AND HIGHER DIMENSIONS 

In this section, we generalize our results to (3 + 1) and 
higher dimensions. We find that there is a one-to-one 
correspondence between (3-t-l) (and higher) dimensional 
string-net condensates and mathematical objects known 
as "symmetric tensor categories." [32] The low energy 
effective theories for these states are gauge theories cou- 
pled to bosonic or fermionic charges. 

Our approach is based on the exactly soluble lattice 
spin Hamiltonian (10). In that model, the spins live on 
the links of the honeycomb lattice. However, the choice 
of lattice was somewhat arbitrary: we could equally well 
have chosen any trivalent lattice in two dimensions. 

Trivalent lattices can also be constructed in three and 
higher dimensions. For example, we can create a space- 
filling trivalent lattice in three dimensions, by "splitting" 
the sites of the cubic lattice (see Fig. 10). Consider the 
spin Hamiltonian (10) for this lattice, where I runs over 
all the vertices of the lattice, and p runs over all the 
"plaquettes" (that is, the closed loops that correspond 
to plaquettes in the original cubic lattice). 

This model is a natural candidate for string-net con- 
densation in three dimensions. Unfortunately, it turns 
out that the Hamiltonian (10) is not exactly soluble on 
this lattice. The magnetic flux operators do not com- 
mute in general. 

This lack of commutativity originates from two differ- 
ences between the plaquettes in the honeycomb lattice 
and in higher dimensional trivalent lattices. The first 
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difference is that in the honeyeomb lattice, neighboring 
plaquettes always share precisely two vertices, while in 
higher dimensions the boundary between plaquettes can 
contain three or more vertices (see Fig. 11). The exis- 
tence of these interior vertices has the following conse- 
quence. Imagine we choose orientation conventions for 
each vertex, so that we have a notion of "left turns" and 
"right turns" for oriented paths on our lattice (such an 
orientation convention can be obtained by projecting the 
3D lattice onto a 2D plane - as in Fig. 11). Then, no 
matter how we assign these orientations the plaquette 
boundaries will always make both left and right turns. 
Thus, we cannot regard the boundaries of the 3D pla- 
quettes as small closed strings the way we did in two 
dimensions (since small closed strings always make all 
left turns, or all right turns). But the magnetic flux op- 
erators Bp only commute if their boundaries are small 
closed strings. It is this inconsistency between the alge- 
braic definition of and the topology of the plaquettes 
that leads to the lack of commutativity. 

To resolve this problem, we need to define a Hamil- 
tonian using the general simple string operators W{dp) 
rather then the small closed strings B^. Suppose 
{colpLvlj), s = 0, 1, ...N are type-s solutions of (18). Af- 
ter picking some "left turn" , "right turn" orientation con- 
vention at each vertex, we can define the corresponding 
type-s simple string operators WsiP) as in (15). Sup- 
pose, in addition, that we choose {u;lj,Ljlj) so that the 
string operators satisfy Wr ■ Wg = J2t ^rstWt (this prop- 
erty ensures that Ws{dp) are analogous to B^). Then, a 
natural higher dimensional generalization of the Hamil- 
tonian (10) is 

N 

H = -J2Qi-J2Wp, Wp = J2asW4dp) (31) 

I p s=0 

For a two dimensional lattice, the conditions (18) are suf- 
ficient to guarantee that the Hamiltonian (31) is an ex- 
actly soluble realization of a doubled topological phase. 
(This is because the plaquette boundaries dp are not 
linked and hence the Ws{dp) all commute). However, 
in higher dimensions, one additional constraint is neces- 
sary. 

This constraint stems from the second, and perhaps 
more fundamental, difference between 2D and higher di- 
mensional lattices. In two dimensions, two closed curves 
always intersect an even mimber of times. For higher 
dimensional lattices, this is not the case. Small closed 
curves, in particular plaquette boundaries, can (in a 
sense) intersect exactly once (see Fig. 11). Because of 
this, the objects must satisfy the additional relation: 

^U=^ij (32) 

One can show that if this additional constraint is satis- 
fied, then (a) the higher dimensional Hamiltonian (31) 
is exactly soluble, and (b) the ground state wave func- 
tion $ is defined by local topological rules analogous to 
(4-7). This means that (31) provides an exactly soluble 
realization of topological phases in (3-1-1) and higher 
dimensions. 




FIG. 11: Tliree plaquettes demonstrating tfic two fundamen- 
tal differences between higher dimensional trivalent lattices 
and the honeycomb lattice. The plaquettes pi, P2 lie in the 
xz plane, while pa is oriented in the xy direction. Notice that 
pi and p2 share three vertices, Ii, I2, la- Also, notice that 
the plaquette boundaries dpi and dpa intersect only at the 
line segment lali- The boundary dpi makes a left turn at 
la, and a right turn at I4. Thus, viewed from far away, these 
two plaquette boundaries intersect exactly once (unlike the 
pair dpi and dp^)- 

Each exactly soluble Hamiltonian is associated with 
a solution {F]^,Z ^^)k^^)k) of (8), (18), (32). By analogy 
with the two dimensional case, we conjecture that there is 
a one-to-one correspondence between topological string- 
net condensed phases in (3 + 1) or higher dimensions, 

and these solutions. The solutions {Fljl^'' ,^jkj^jk) '-or- 
respond to a special class of tensor categories - symmet- 
ric tensor categories. [32] Thus, just as tensor categories 
are the mathematical objects underlying string conden- 
sation in (2 + 1) dimensions, symmetric tensor categories 
are fundamental to string condensation in higher dimen- 
sions. 

There arc relatively few solutions to (8), (18), (32). 
Physically, this is a consequence of the restrictions on 
quasiparticle statistics in 3 or higher dimensions. Un- 
like in two dimensions, higher dimensional quasiparticles 
necessarily have trivial mutual statistics, and must be ei- 
ther bosonic or fermionic. From a more mathematical 
point of view, the scarcity of solutions is a result of the 
symmetry condition (32). Doubled topological phases, 
such as Chern-Simons theories, typically fail to satisfy 
this condition. 

However, gauge theories do satisfy the symmetry con- 
dition (32) and therefore do correspond to higher dimen- 
sional string- net condensates. Recall that the gauge the- 
ory solution to (8) is obtained by (a) letting the string- 
type index i run over the irreducible representations of 
the gauge group, (b) letting the numbers di be the dimen- 
sions of the representations, and (c) letting the 6 index 
object Fjfi^ be the 6j symbol of the group. One can 
check that this also provides a solution to (18), (32), if 
we set ijjli.^^^^ = — 1 when j = k and the invariant ten- 
sor in fc ® fc ® i* is antisymmetric in the first two indices, 
and J. = 1 otherwise. This result is to be expected, 
since the string-net picture of gauge theory (section II) 
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is valid in any number of dimensions. Thus, it is not 
surprising that gauge theories can emerge from higher 
dimensional string-net condensation. 

There is another class of higher dimensional string- 
net condensed phases that is more interesting. The low 
energy effective theories for these phases are variants of 
gauge theories. Mathematically, they are obtained by 
twisting the usual gauge theory solution by 

(i},=a;},.(-l)^(^)^W (33) 

Here P{i) is some assignment of parity ( "even" , or "odd" ) 
to each representation i. The assignment must be self- 
consistent in the sense that the tensor product of two 
representations with the same (different) parity, decom- 
poses into purely even (odd) representations. If all the 
representations are assigned an even parity - that is P{i) 
is trivial - then the twisted gauge theory reduces to stan- 
dard gauge theory. 

The major physical distinction between twisted gauge 
theories and standard gauge theories is the quasiparti- 
cle spectrum. In standard gauge theory, the fundamen- 
tal quasiparticles are the electric charges created by the 
A'' -|- 1 string operators Wj. These quasiparticles are all 
bosonic. In contrast, in twisted gauge theories, all the 
quasiparticles corresponding to "odd" representations i 
are fermionic. 

In this way, higher dimensional string-net condensa- 
tion naturally gives rise to both emerging gauge bosons 
and emerging fermions. This feature suggests that gauge 
interactions and Fermi statistics may be intimately con- 
nected. The string- net picture may be the bridge between 
these two seemingly unrelated phenomena. [29] 

In fact, it appears that gauge theories coupled to 
fermionic or bosonic charged particles are the only pos- 
sibilities for higher dimensional string-net condensates: 
mathematical work on symmetric tensor categories sug- 
gests that the only solutions to (8), (18), (32) are those 
corresponding to gauge theories and twisted gauge theo- 
ries. [48] 

We would like to point out that (3 -|- 1) dimensional 
string-net condensed states also exhibit membrane con- 
densation. These membrane operators are entirely anal- 
ogous to the string operators. Just as open string oper- 
ators create charges at their two ends, open membrane 
operators create magnetic flux loops along their bound- 
aries. Furthermore, just as string condensation makes the 
string unobservable, membrane condensation leads to the 
unobservability of the membrane. Only the boundary of 
the membrane - the magnetic flux loop - is observable. 

VI. EXAMPLES 
A. N = 1 string model 

We begin with the simplest string-net model. In the 
notation from section III, this model is given by 

1. Number of string types: N = 1 

2. Branching rules: (no branching) 



3. String orientations: 1* = 1. 

In other words, the string-nets in this model contain one 
unoriented string type and have no branching. Thus, 
they are simply closed loops. (See Fig. 3a). 

We would like to find the different topological phases 
that can emerge from these closed loops. According to 
the discussion in section IV, each phase is captured by 
a fixed-point wave function, and each fixed-point wave 
function is specified by local rules (4-7) that satisfy the 
self-consistency relations (8). It turns out that (8) have 
only two solutions in this case (up to rescaling): 

do = 1 

di = FllS = ±1 

pOOO _ plOl _ pOll _ 1 
^000 — ^101 — ^011 — ^ 

pooo _ piio _ pioi _ jpon _ 1 fyA-^ 
— ^001 — ^010 — ^100 — ^ \'^^) 

where the other elements of F all vanish. The corre- 
sponding local rules (4-7) are: 




*(■><■)= ± ^ (■ I) (35) 

We have omitted those rules that can be derived from 
topological invariance (4). 

The fixed-point wave functions $± satisfying these 
rules are given by 

$±(X) = (±1)^= (36) 

where Xc is the number of disconnected components in 
the string configuration X. 

The two fixed-point wave functions $± correspond 
to two simple topological phases. As we will sec, 
corresponds to Z2 gauge theory, while is a U{1) x 
^(1) Chern-Simons theory. (Actually, other topological 
phases can emerge from closed loops - such as in Ref. [42- 
44]. However, we regard these phases as emerging from 
more complicated string-nets. The closed loops organize 
into these effective string-nets in the infrared limit). 

The exactly soluble models (10) realizing these two 
phases can be written as spin 1/2 systems with one spin 
on each link of the honeycomb lattice (see Fig. 12). We 
regard a link with cr^ = — 1 as being occupied by a type-1 
string, and the state cr^ = +1 as being unoccupied (or 
equivalently, occupied by a type-0 or null string). The 
Hamiltonians for the two phases are of the form 

I p 

The electric charge term is the same for both phases 
(since it only depends on the branching rules): 

Q/,±=^(i+ n (37) 

legs of / 
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FIG. 12: The Hamiltonians (39), (40), realizing the two iV = 1 
string-condensed phases. Each circle denotes a spin-1/2 spin. 
The links with — —1 are thought of as being occupied by 
a type-1 string, while the links with cr^ = +1 are regarded as 
empty. The electric charge term acts on the three legs of the 
vertex / with . The magnetic energy term acts on the 6 
edges of the plaquette p with ct^ , and acts on the 6 legs of p 
with an operator of the form /(cr^). For the Z2 phase, / = 1, 
while for the Chern-Simons phase, f{x) 



,-(l-^)/2 



The magnetic terms for the two phases are 



1 



Bp.± - :j(Bp ±±Bp^) 



n ' 

edges of p 




n (^/^)- 

legs of p 



(38) 



Pr, 



The 



where Pp is the projection operator Pp = Yiiep 
projection operator Pp can be omitted without affecting 
the physics (or the exact solubility of the Hamiltonian). 
We have included it only to be consistent with (10). If 
we omit this term, the Hamiltonian for the first phase 
(<&+) reduces to the usual exactly soluble Hamiltonian of 
Z2 lattice gauge theory (neglecting numerical factors): 



E n E n ' 

I legs of 7 p edges of p 

The Hamiltonian for the second phase, 



(39) 



^---E n -f+E( n n 

I legs of J p edges of p legs of p 

(40) 

is less familiar. However, one can check that in both 
cases, the Hamiltonians are exactly soluble and the two 
ground state wave functions are precisely $± (in the 
basis) . 

Next we find the quasiparticle excitations for the two 
phases, and the corresponding S-matrix and twists Oa- 

In both cases, equation (22) has 4 irreducible solu- 
tions {riai.sT^^a stT^^a st)j = l!2,3,4 - Corresponding to 
4 quasiparticle types. For the first phase {^+) these so- 
lutions are given by: 



ni,o = 1, "1,1 = 0, 

n2fi = 0, n2,i = 1, 

713,0 = 1, "3,1 = 0, 

714,0 = 0, 774,1 = 1, 



"1,000 
"2.110 
"3.000 
"4,110 



= Ij ^1,001 = 1 

— ^2,111 = 1 

— 1' ^3,001 = ^1 
= 1, in = -1 



FIG. 13: A closed string operator W{P) for the two models 
(39), (40). The path P is drawn with a thick line, while the legs 
are drawn with thin lines. The action of the string operators 
(41), (44) on the legs is different for legs that branch to the 
right of P, "R-legs" , and legs that branch to the left of P, "L- 
legs." Similarly, we distinguish between "R-vertices" and "L- 
vertices" which are ends of "R-leg" and "L-leg" respectively. 



The other elements of f2 vanish. In all cases = il. 
The corresponding string operators for a path P are 



Id 



n 



edges of P 



R-lcgs 

n -I n -I 

edges of P R-legs 



(41) 



where the "R-legs" k are the legs that are to the right 
of P. (See Fig. 13). Technically, we should multiply 
these string operators by an additional projection opera- 
tor Wj^pQi, in order to be consistent with the general 
result (19). However, we will neglect this factor since it 
doesn't affect the physics. 

Once we have the string operators, we can easily cal- 
culate the twists and the S-matrix. We find: 



(42) 
(43) 



This is in agreement with the twists and S-matrix for 
Z2 gauge theory: Wi creates trivial quasiparticles, W2 
creates magnetic fluxes, Wz creates electric charges, W4 
creates electric/magnetic bound states. 
In the second phase ($_), we flnd 




"■1,0 


= 1, 


"1,1 


= 0, 


0" 

"1.000 


= 1, 


"1,001 ~ 


»^2,0 


= 0, 


«2,1 


= 1, 


"2.110 


= 1, 


^2,111 ~ * 


'^3,0 


= 0, 


«3,1 


= 1, 


"3,110 


= 1, 


^3,111 — 


»^4,0 


= 1, 


«4,1 


= 0, 


"4,000 


= 1, 


"4,001 ~ ^ 



Once again, the other elements of Vl vanish. Also, in all 
cases, Ct = W . The corresponding string operators for a 
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path P are 
Wi = Id 

^2 = n n n i-^r 

edges of P R-legs L-vertices 

^3 = n -J n (-^)^ n (-^r 

edges of P R-legs 



L-vertices 



W4 



n 

R-legs 



(44) 



where the "L-vertices" I are the vertices of P adjacent to 
legs that are to the left of P. The exponent s/ is defined 
by Si = j{l — cr?)(l + cr|), where i, j are the links just 
before and just after the vertex I, along the path P. (See 
Fig. 13). 

We find the twists and S-matrix are 



1, 



1 



5 = 




(45) 
(46) 



We see that Wi creates trivial quasiparticles, W2, W3 
create semions with opposite chiralities and trivial mu- 
tual statistics, and W4 creates bosonic bound states of 
the semions. These results agree with the f/(l) x U{1) 
Chern-Simons theory 



L = —Kijai^d^ajxe'^''^, 



I, J =1,2 



with if -matrix 



K 



2 
-2 



(47) 



(48) 



Thus the above C/(l) x U{1) Chern-Simons theory is the 
low energy effective theory of the second exactly soluble 
model (with di = —1). 

Note that the Z2 gauge theory from the first exactly 
soluble model (with di = 1) can also be viewed as a 
C/(l) X U{1) Chern-Simons theory with if-matrix [23] 



K = 



2 
2 



(49) 



B. N = 1 string- net model 

The next simplest string-net model also contains only 
one oriented string type - but with branching. Simple 
as it is, we will see that this model contains non-Abclian 
anyons and is theoretically capable of universal fault tol- 
erant quantum computation [33] . Formally, the model is 
defined by 

1. Number of string types: N = 1 

2. Branching rules: {{1, 1, 1}} 

3. String orientations: 1* = 1. 



The string- nets are unoricntcd trivalcnt graphs. To find 
the topological phases that can emerge from these ob- 
jects, we solve the self-consistency relations (8). We find 
two sets of self-consistent rules: 



O . $ 



=7± 



1/2 



m> <m 



(50) 



where 7± = ■ (Once again, we have omitted those 
rules that can be derived from topological invariance). 
Unlike the previous case, there is no closed form expres- 
sion for the wave function amplitude. 

Note that the second solution, di = does not 

satisfy the unitarity condition (14). Thus, only the first 
solution corresponds to a physical topological phase. As 
we will see, this phase is described by an 503(3) x 503(3) 
Chern-Simons theory. 

As before, the exactly soluble realization of this phase 
(10) is a spin-1/2 model with spins on the links of the 
honeycomb lattice. We regard a link with cr^ = —1 as 
being occupied by a typc-1 string, and a link with = 1 
as being unoccupied (or equivalcntly occupied by a type- 
string). However, in this case we will not explicitly 
rewrite (10) in terms of Pauli matrices, since the resulting 
expression is quite complicated. 

We now find the quasiparticles. These correspond to 
irreducible solutions of (22). For this model, there are 4 
such solutions, corresponding to 4 quasiparticles: 



1 : ni,o = 1, ni,i = 0, 000 = ^1, 

2 : ^2,0 = 0, 712,1 = 1, ^ 



001 



1 



2.110 



(51) 



7+ e , "2,111 ^7-1- e 



3 : ^3,0 = 0, 713,1 = 1, il3,iio = 1, 



"3,111 

4 : 774,0 — 1, 'T-4,1 — ^, ^'4 000 



"4.001 



7+ I "3,111 

^' "4,000 ~ ^' ^'4,110 

7+ 



^-l/2^_3W5 



2 qO Ql -^-5/2 

"4,111—7+ ; "4,111—7+ ) 



f^iioi = (^^ioii)* = 7;"/'(2 ■ 



„37ri/5 



■7+e 



-37ri/5 



In all cases, = 

We can calculate the twists and the S-matrix. We find: 

g«i ^ gi02 = e-*'''/^, = e^'^'/^, e'^* = 1 (52) 



(53) 



We conclude that Wi creates trivial quasiparticles, W2, 
W3 create (non-Abelian) anyons with opposite chirali- 
ties, and W4 creates bosonic bound states of the anyons. 
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These results agree with 6*03 (3) x SO3 (3) Chern-Simons 
theory, the so-cahed doubled "Yang-Lee" theory. 

Researchers in the field of quantum computing have 
shown that the Yang- Lee theory can function as a univer- 
sal quantum computer - via manipulation of non-Abelian 
anyons. [33] Therefore, the spin-1/2 Hamiltonian (10) as- 
sociated with (50) is a theoretical realization of a univer- 
sal quantum computer. While this Hamiltonian may be 
too complicated to be realized experimentally, the string- 
net picture suggests that this problem can be overcome. 
Indeed, the string-net picture suggests that generic spin 
Hamiltonians with a trivalent graph structure will ex- 
hibit a Yang-Lee phase. Thus, much simpler spin-1/2 
Hamiltonians may be capable of universal fault tolerant 
quantum computation. 



C. N = 2 string-net models 

In this section, we discuss two N = 2 string-net mod- 
els. The first model contains one oriented string and its 
dual. In the notation from section III, it is given by 

1. Number of string types: N = 2 

2. Branching rules: {{1, 1, 1}, {2, 2, 2}} 

3. String orientations: 1* = 2, 2* = 1. 

The string-nets are therefore oriented trivalent graphs 
with Z3 branching rules. The string-net condensed 

phases correspond to solutions of (8) . Solving these equa- 
tions, we find two sets of self-consistent local rules: 

$ = ± $ 

^ (■><■) =^/^^•^(■I■) (54) 

The corresponding fixed-point wave functions $± are 
given by 



$±(x) = (V±T)'''<=-^"/' 



(55) 



where Xc, Xy, are the number of connected components, 
and vertices, respectively in the string- net configuration 
X. As before, we can construct an exactly soluble Hamil- 
tonians, find the quasiparticles for the two theories and 
compute the twists and S-matrices. We find that the first 
theory is described by a gauge theory, while the 
second theory $_ is described by a U{1) x U{1) Chern- 
Simons theory with JT-matrix 



K 



3 
-3 



Both theories have 3^ = 9 elementary quasiparticles. 

In the case of Z3, these quasiparticles are electric 
charge/magnetic flux bound states formed from the 3 
types of electric charges and 3 types of magnetic fluxes. 



In the case of the Chern-Simons theory, the quasiparti- 
cles are bound states of the two fundamental anyons with 
statistical angles ±7r/3. 

The flnal example we will discuss contains two unori- 
ented strings. Formally it is given by 

1. Number of string types: N = 2 

2. Branching rules: {{1, 2, 2}, {2, 2, 2}} 

3. String orientations: 1* = 1, 2* = 2. 

The string-nets are unoriented trivalent graphs, with 
edges labeled with 1 or 2. We find that there is only 
one set of self-consistent local rules: 



--di ■ $ 



■nil' 



2 i 2 
1 1 2 



n=0 



22m . ^ 

22n ^ 



mm 



(56) 



where do = rfi = 1, ^2 = 2, and F'^™ is the matrix 




77.22m 



If we construct the Hamiltonian (10), we find that it 
is equivalent to the standard exactly lattice gauge theory 
Hamiltonian [20] with gauge group 5*3 - the permuta- 
tion group on 3 objects. One can show that this theory 
contains 8 elementary quasiparticles (corresponding to 
the 8 irreducible representations of the quantum double 
D{S-i)). These quasiparticles are combinations of the 3 
electric charges and 3 magnetic fiuxes. 



VII. CONCLUSION 

In this paper, we have shown that quantum systems 
of extended objects naturally give rise to topological 
phases. These phases occur when the extended objects 
(e.g. string-nets) become highly fluctuating and con- 
dense. This physical picture provides a natural mech- 
anism for the emergence of parity invariant topological 
phases. Microscopic degrees of freedoms (such as spins 
or dimers) can organize into effective extended objects 
which can then condense. We hope that this physical 
picture may help direct the search for topological phases 
in real condensed matter systems. It would be interest- 
ing to develop an analogous picture for chiral topological 
phases. 
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We have also found the fundamental mathematical 
framework for topological phases. We have shown that 
each (2 + 1) dimensional doubled topological phase is as- 
sociated with a 6 index object F^j^ and a set of real 
numbers di satisfying the algebraic relations (8). All 
the universal properties of the topological phase are con- 
tained in these mathematical objects (known as tensor 
categories). In particular, the tensor category directly 
determines the quasiparticle statistics of the associated 
topological phase (28, 30). This mathematical frame- 
work may also have applications to phase transitions and 
critical phenomena. Tensor categories may characterize 
transitions between topological phases just as symmetry 
groups characterize transitions between ordered phases. 

We have constructed exactly soluble (2 + 1)D lattice 
spin Hamiltonians (10) realizing each of these doubled 
topological phases. These models unify (2 + 1)D lattice 
gauge theory and doubled Chern-Simons theory. One 
particular Hamiltonian - a realization of the doubled 
Yang-Lcc theory - is a spin 1/2 model capable of fault 
tolerant quantum computation. 

In higher dimensions, string-nets can also give rise to 
topological phases. However, the physical and mathe- 
matical structure of these phases is more restricted. On 
a mathematical level, each higher dimensional string-net 
condensate is associated with a special kind of tensor 
category - a symmetric tensor category (18), (32). More 
physically, we have shown that higher dimensional string- 
net condensation naturally gives rise to both gauge inter- 
actions and Fermi statistics. Viewed from this perspec- 
tive, string-net condensation provides a mechanism for 
unifying gauge interactions and Fermi statistics. It may 
have applications to high energy physics [30] . 

From a more general point of view, all of the phases 
described by Landau's symmetry breaking theory can 
be understood in terms of particle condensation. These 
phases are classified using group theory and lead to emer- 
gent gapless scalar bosons [49, 50], such as phonons, spin 
waves, etc . In this paper, we have shown that there is a 
much richer class of phase - arising from the condensation 
of extended objects. These phases are classified using 
tensor category theory and lead to emergence of Fermi 
statistics and gauge excitations. Clearly, there is whole 
new world beyond the paradigm of symmetry breaking 
and long range order. It is a virgin land waiting to be 
explored. 

We would like to thank Pavel Etingof and Michael 
Freedman for useful discussions of the mathematical as- 
pects of topological field theory. This research is sup- 
ported by NSF Grant No. DMR-01-23156 and by NSF- 
MRSEC Grant No. DMR-02-13282. 



APPENDIX A: 



GENERAL STRING-NET 
MODELS 



In this section, we discuss the most general string-net 

models. These models can describe all doubled topo- 
logical phases, including all discrete gauge theories and 
doubled Chern-Simons theories. 



In these models, there is a "spin" degree of freedom at 
each branch point or node of a string-net, in addition to 
the usual string-net degrees of freedom. The dimension 
of this "spin" Hilbert space depends on the string types 
of the 3 strings incident on the node. 

To specify a particular model one needs to provide 

a 3 index tensor Sijk which gives the dimension of the 
spin Hilbert space associated with k} (in addition to 
the usual information). The string-net models discussed 
above correspond to the special case where 6ijk = 0, 1 for 
all k. In the case of gauge theory, 6ijk is the number 
of copies of the trivial representation that appear in the 
tensor product i®j®k. Thus we need the more general 
string-net picture to describe gauge theories where the 
trivial representation appears multiple times in i® j^k. 

The Hilbert space of the string-net model is defined 
in the natural way: the states in the string-net Hilbert 
space are linear superpositions of different spatial config- 
urations of string-nets with different spin states at the 
nodes. 

One can analyze string-net condensed phases as before. 
The universal properties of each phase are captured by 
a fixed-point ground state wave function The wave 
function $ is specified by the local rules (4), (5) and 
simple modifications of (6), (7): 




The complex numerical constant Fj^j^ is now a complex 
tensor {Flj^)'^l of dimension Sijm x Skim* x <^in( x Sjkn' • 
One can proceed as before, with self-consistency con- 
ditions, fixed-point Hamiltonians, string operators, and 
the generalization to (3-1-1) dimensions. The exactly sol- 
uble models are similar to (10). The main difference is 
the existence of an additional spin degree of freedom at 
each site of the honeycomb lattice. These spins account 
for the degrees of freedom at the nodes of the string-nets. 



APPENDIX B: SELF-CONSISTENCY 
CONDITIONS 



In this section, we derive the self-consistency condi- 
tions (8). We begin with the last relation, the so-called 
"pentagon identity", since it is the most fundamental. 
To derive this condition, we use the fusion rule (7) 



jk 



to relate the amplitude $ \3^^i^/ ] to the amplitude 

V A fm J 

^ in two distinct ways (see Fig. 14). On the 

one hand, we can apply the fusion rule (7) twice to obtain 
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FIG. 14: The fusion rule (7) can be used to relate the am- 
plitude of (a) to the amplitude of (c) in two different ways. 
On the one hand, we can apply the fusion rule (7) twice - 
along the links denoted by solid arrows - to relate (a) —> (b) 
(c). But we can also apply (7) three times - along the links 
denoted by dashed arrows - to relate (a) (d) (e) (c). 
Self-consistency requires that the two sequences of the opera- 
tion lead to the same linear relations between the amplitudes 
of (a) and (c). 




r 

Epjip prif ^ 
q*kr* mis* 




(Here, we neglected to draw a shaded region surround- 
ing the whole diagram. Just as in the local rules (4-7) 
the ends of the strings i, j, k, I, m arc connected to some 
arbitrary string- net configuration). But we can also ap- 
ply the fusion rule (7) three times to obtain a different 
relation: 



n 

n,s \ 1/ 'm 

n.r.K \ i //m J 



If the rules are self-consistent, then these two relations 
must agree with each other. Thus, the two coefficients 





of $ 




must be the same. This equality implies 



the pentagon identity (8). 

The first two relations in (8) are less fundamental. In 
fact, the first relation is not required by self- consistency 
at all; it is simply a useful convention. To see this, con- 
sider the following rcscaling transformation on wave func- 
tions $ ^ Given a string- net wave function $, we can 
obtain a new wave function $ by multiplying the ampli- 
tude ^{X) for a string-net configuration X by an arbi- 
trary factor f{i,j,k) for each vertex {i,j,k} in X. As 
long as f{i, j, k) is symmetric in i,j,k and /(O, = 1, 




-m J 



k 

(d) 



FIG. 15: Four string-net configurations related by tetrahe- 
dral symmetry. In diagram (a), we show the tetrahedron cor- 
responding to Gjf,^. In diagrams (b), (c),(d), we show the 

tetrahedrons Gf 



'J^n , GjlJ^,, G^J^^i, obtained by reflecting (a) 
in 3 different planes: the plane joining n to the center of m, 
the plane joining m to the center of n, and the plane joining 
i to the center of k. The four tetrahedrons correspond to the 
four terms in the second relation of (8). 



this operation preserves the topological invariance of $. 
The rescaled wave function $ satisfies the same set of 



local rules with rescaled F, 



^kln 



fpijm 



Tpljm 
^kln 



f{hj,m)f{k,l,m*) 
f{n,l,i)f{j,k,n*) 



(Bl) 



Since $ and i> describe the same quantum phase, we 
regard F and F as equivalent local rules. Thus the first 
relation in (8) is simply a normalization convention for F 
or $ (except when i,j or k vanishes; these cases require 
an argument similar to the derivation of the pentagon 
identity) . 

The second relation in (8) has more content. This re- 
lation can be derived by computing the amplitude for a 
tetrahedral string-net configuration. We have: 



-^kln ^fc'I' 





=F, 



ijm TjnkT Tpn'i'l* 



dkdidn^{$) 



(B2) 
(B3) 



We define the above combination in the front of $(0) 



as: 



Pkln ViVjVkVl 



(B4) 



Imagine that the above string-net configuration lies on 
a sphere. In that case, topological invariance (together 
with parity invariance) requires that G^^^ be invariant 
under all 24 symmetries of a regular tetrahedron. The 
second relation in (8) is simply a statement of this tetra- 
hedral symmetry requirement - written in terms of -F^;™- 
(See Fig. 15). 

In this section, we have shown that the relations (8) 
arc necessary for self-consistency. It turns out that these 
relations are also sufficient. One way of proving this is to 
use the lattice model (10). A straightforward algebraic 
calculation shows that the ground state of 10 obeys the 
local rules (4-7), as long as (8) is satisfied. This estab- 
lishes that the local rules are self-consistent. 
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FIG. 16: The fattened honeycomb lattice. The strings are 
forbidden in the shaded region. A string state in the fattened 
honeycomb lattice (a) can be viewed as a superposition of 
string states on the links (b). 



APPENDIX C: GRAPHICAL REPRESENTATION 
OF THE HAMILTONIAN 

In this section, wc provide an alternative, graphical, 
representation of the lattice model (10). This graphical 
representation provides a simple visual technique for un- 
derstanding properties (a)-(c) of the Hamiltonian (10). 

We begin with the 2D honeycomb lattice. Imagine we 
fatten the links of the lattice into stripes of finite width 
(see Fig. 16). Then, any string- net state in the fat- 
tened honeycomb lattice (Fig. 16a) can be viewed as a 
superposition of string-net states in the original, imfat- 
tened lattice (Fig. 16b). This mapping is obtained via 
the local rules (4-7). Using these rules, we can relate the 
amplitude $(X) for a string- net in the fattened lattice to 
a linear combination of string-net amplitudes in the orig- 
inal lattice: ^{X) = ^ai$(Xi). This provides a natural 
linear relation between the states in the fattened lattice 



and those in the imfattcncd lattice: \X) = Y^ai\Xi). 
This linear relation is independent of the particular way 
in which the local rules (4-7) are applied, as long as the 
rules are self-consistent. 

In this way, the fattened honeycomb lattice provides 
an alternative notation for representing the states in the 
Hilbert space of (10). This notation is useful because 
the magnetic energy operators are simple in this rep- 
resentation. Indeed, the action of the operator on 



" ) is equivalent to simply 



the string-net state 

'\ /i 
f 

adding a loop of type-s string: 





As we described above, we can use the local rules (4-7) 
to rewrite ^•Ai'' Yz-o ) as a linear combination of the 




physical string- net states with strings only on the links, 
that is to reduce Fig. 17a to Fig. 17b. This allows us to 
obtain the matrix elements of B^. 

The following is a particular way to implement the 
above procedure: 




Notice that (CI) is exactly (13). Thus, the graphical 
representation of Bp agrees with the original algebraic 

definition. 

Using the graphical representation of B^ we can eas- 
ily show that Bp\ and commute. The derivation 



is much simpler then the more straightforward algebraic 
calculation. First note that these operators will commute 
if pi, P2 are well-separated. Thus, we only have to con- 
sider the case where pi and p2 are adjacent, or the case 
where pi, p2 coincide. We begin with the nearest neigh- 
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FIG. 17: The action of Bp is equivalent to adding a loop of 
type-s string. The resulting string-net state (a) is actually 
a linear combination of the string-net states (b). The coef- 
ficients in this linear relation can be obtained by using the 
local rules (4-7) to reduce (a) to (b). 
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FIG. 18: The action of Bp\B!pi^ on the string-net state (a) 
can be represented by adding two loops of type-si and type- 
S2 strings as shown in (b). The string-net state (b) is a linear 
combination of the string-net states (c). The coefficients are 
obtained by using (4 -7) to reduce (b) to (c). 



bor case. The action of B^^B^^ on the string- net state 
Fig. 18a can be represented as Fig. 18b. Fig. 18b can 

then be related to a hncar combination of the string-net 
states shown in Fig. 18c. The coefficients in this relation 
are the matrix elements of Bt^ Bt? . But by the same ar- 
gument, the action of B^^B^^^ can also be represented by 
Fig. 18b. We conclude that Bt? Bt} , Bt} Bt? have the 

^ P2 Pl ' Pi P2 

same matrix elements. Thus, the two operators commute 
in this case. 

On the other hand, when = P2, we have 



b;'b;^ 




— Z^-^s*s2fe*-^s;feo '*« 

k 




*S2Sl 




Thus, 



p 



y^Jk'S2SiB^. 



(C2) 



(C3) 



Since (Jfc*s2Si is symmetric in S2, Si, we conclude that 
Bp^ B^ = B^B^ , so the operators commute in this case 
as well. This establishes property (a) of the Hamiltonian 
(10). 

Equation (C3) also sheds light on the spectrum of the 
Bp operators. Let the simultaneous eigenvalues of Bp 
(with p fixed) be {bg}. Then, by (C3) these eigenvalues 
satisfy 

k 

We can view this as an eigenvalue equation for the {N + 
l)x{N + 1) matrix Ms^, defined by Mj^_^. = 5j*s2i- The 



simultaneous eigenvalues are simply the simultaneous 
eigenvalues of the matrices M^^ . In particular, this means 
that the index q ranges over a set of size A'' -I- 1. 

Each value of q corresponds to a different possible state 
for the plaquette p. The magnetic energies of these -I- 1 
different states are given by: Eg = — J2s ^sbg- Depend- 
ing on the parameter choice as, all on the plaquettes p 
will be in one of these states q. In this way, the Hamilto- 
nian (10) can be in A''-!- 1 different quantum phase. This 
establishes property (b) of the Hamiltonian (10). 

One particular state q is particularly interesting. This 
state corresponds to the simultaneous eigenvalues = 
ds- It is not hard to show that the parameter choice 
makes this state energetically favorable. In 



fact, using (03) one can show that Bp is a projector for 
this parameter choice, and that Bp = 1 for this state. 

Furthermore, the ground state wave function for this 
parameter choice obeys the local rules (4-7). One way 



to see this is to compare Bp 
For the first state, we find 



with Br 



Br 



s 

E 



s 



3-s 



Ekdl 



^ s*sj* 



s 



-E 



VjVs 



s 



For the second state, we find the same result: 



Br 



It follows that 



E 



ViVs 



j 




so 



This result means that the strings can be moved through 
the forbidden regions at the center of the hexagons. 
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FIG. 19: The action of the string operator Wa{P) is equiva- 
lent to adding a type-a string along the path P. The result- 
ing string-net state can be reduced to a Uucar combination of 
states on the honeycomb lattice, using the local rules (4-7), 
(Dl). 



Thus, the local rules which were originally restricted to 
the fattened honeycomb lattice can be extended through- 
out the entire 2D plane. The wave function $ obeys these 
continuum local rules and has a smooth continuum limit. 
We call such a state smooth topological state. This es- 
tablishes property (c) of the Hamiltonian (10). 

The wave fimctions of some smooth topological states 
are positive definite. So those wave functions can be 
viewed as the statistical weights of certain statistical 
models in the same spatial dimensions. What is inter- 
esting is that those statistical models are local models 
with short-ranged interactions [10, 35, 36]. 



APPENDIX D: GRAPHICAL REPRESENTATION 
OF THE STRING OPERATORS 

In this section, we describe a graphical representation 
of the long string operators Wa{P)- Just as in the pre- 
vious section, this representation involves the fattened 
honeycomb lattice. The action of the string operator 
Wa{P) on a general string state X, is simply to create a 
string labeled a along the path P. The resulting string- 
net state can then be reduced to a linear combination of 
string-net states on the unfattened lattice. The coeffi- 
cients in this linear combination are the matrix elements 
oiWo^iP). 

However, none of the rules (4-7) involve strings labeled 



a, nor do they allow for crossings. Thus, the reduction to 
string-net states on the unfattened lattice requires new 
local rules. These new local rules are defined by the 4 
index objects j-, and the integers na,i- 



i 



3 St 



jst 



(Dl) 



Here, cr, r are the two indices of the matrix Q^j. (Until 
now, we've neglected to write out these indices explic- 
itly). 

After applying these rules, we then need to join to- 
gether the resulting string-nets. The "joining rule" for 
two string types si, S2 is as follows. If si ^ S2, we don't 
join the two strings: we simply throw away the diagram. 
If Si = S2, then we join the two strings and contract the 
two corresponding indices cti, (T2- That is, we multiply 
the two matrices together in the usual way. Using the 
same approach as (CI), one can show that the graphical 
definition of Wa{P) agrees with the algebraic definition 
(19). 

In the previous section, we used the graphical repre- 
sentation of Bp to show that these operators commute. 
The string operators Wa {P) can be analyzed in the same 
way. With a simple graphical argument one can show 
that the string operators Wa (P) commute with the mag- 
netic operators provided that (4-7), (Dl) satisfy the 
conditions 



j>7 a 



(D2) 
(D3) 



These relations are precisely the commutativity condi- 
tions (22), written in graphical form. 
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